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An alternative approach to lattice gauge theory has been under development for the past decade. It is based on 
discretizing the operator Heisenberg equations of motion in such a way as to preserve the canonical commutation 
relations at each lattice site. It is now known how to formulate a non-Abelian gauge theory within this framework. 
The formulation appears to be free of fermion doubling. Since the theory is unitary, a time-development operator 
(Hamiltonian) can be constructed. 



1. QUANTUM MECHANICS 

Let us consider a one-dimensional field theory, 
namely quantum mechanics, governed by a con- 
tinuum Hamiltonian H — p 2 /2 + V(q). The cor- 
responding Heisenberg equations are 



V = -V'(q), 



q = p. 



(1) 



These equations are to be solved subject to the 
initial condition [q(0),p(0)] = i. The equations 
(|j) preserve unitarity in the sense that at any 
subsequent time t the canonical commutation re- 
lations continue to hold, [q(t),p(t)] = i. 

1.1. Finite-element discretization 

We discretize the above equations of motion 
by dividing the time interval [0, T] into N equal 
intervals of length h. On each interval we define 
a local time variable t, < t < h, and write the 
dynamical variables as rth degree polynomials in 
t, 



(2) 



k=0 



k=0 



We determine the 2r + 2 operator coefficients, etfc, 
6^, by imposing continuity between the intervals, 
that is, at the lattice sites, and by imposing the 
equations (Q) at r different points in the interval, 
at U = aih. Unitarity, that is, that = i 

exactly at each lattice site, then uniquely deter- 
mines that the points at which the equations of 
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motion are to be imposed are the Gaussian knots, 
the solutions of P r (2a — 1) = 0. Remarkably, 
these are precisely the points at which classically 
the numerical error is minimized, the relative er- 
ror after TV steps then going like N~ 2r . 

1.2. Time-evolution operator 

Because the lattice theory is unitary, there 
must exist a unitary time-evolution operator 
which advances operators through the lattice, 



q n+1 = Uq n U~ 



Pn 



U Pn U~ 



(3) 



For linear finite elements, r — 1, U is given by 



U 



e ihp 2 n /i e ihA(q n ) e ihpl I A 



(4) 



A{q) = ^[q-g-^Aq/h^ + V^iAq/h 2 )]^) 



g{x) = —x + V\x). 



(6) 



It is apparent that this is an implicit scheme, as 
the equations of motion must be solved to con- 
struct the operator A. Of course, we can use 
this evolution operator to define a self-adjoint 
lattice Hamiltonian, H = — j-lnU, which differs 
from the continuum Hamiltonian by 0(h 2 ). [The 
order r finite element prescription would give a 
Hamiltonian differing from the continuum one by 
0(h 2r )}. The details are given in @. 

1.3. Harmonic-Oscillator Matrix Elements 
of the Time-Evolution Operator 

Remarkably, in view of the above, a closed ex- 
pression, involving g not <7 _1 , can be derived for 



2 



matrix elements of the transfer operator in a har- 
monic oscillator basis Although this for- 
mula is fairly complex, it can be easily used to 
extract expansions in powers of the lattice spacing 
h. For the anharmonic oscillator with potential 
V = Xq 2k /2k, we hnd in the harmonic oscillator 
groundstate 

(U) = l + ihX^+Vfia) - ^A 2 /( fe+1 ) S (a), (7) 

where a = Xj 2k+2 , 7 is the width of the harmonic 
oscillator state, and 



/(«) 



s(a 



1 



4a 1 /('=+i) 
1 



3 -4a 



2k- lT(k+ 1/2) 



16a 2/(fe+l 

4a 2 T(2fc + 1/2) 



r(i/2) 



+ 



k 2 r(i/2) 



(9) 



We match this expansion to that of exp(iuioh) 
to determine approximations to the ground state 
energy. If we use the order h data only, we must 
determine the variational parameter a by, say, 
the principle of minimum sensitivity. The energy 
so determined is accurate to a few percent. If 
the order h 2 data is used a is determined, and 
the results are considerably more accurate. If, 
instead, we use a two-state approximation to the 
anharmonic oscillator ground state, improvement 
by more than an order of magnitude results. See 
Table [l]. Wavefunctions are equally impressive. 

2. NON-ABELIAN GAUGE THEORY 

Recently it was discovered how to implement 
non-Abelian gauge invariance on a linear finite- 
element lattice y. This was done construc- 
tively, by gauging the free lattice Dirac and 
Yang-Mills equations (which, by the linear finite- 
element prescription, are forward difference equa- 
tions, forward-averaged in the non-differentiated 
directions), and simultaneously determining the 
form of the gauge transformation. The result is 
most conveniently expressed in terms of the fol- 
lowing link operator defined on a particular finite- 
element (hypercube): 



where m is a local four-vector coordinate, having 
components or 1, indicating the 16 corners of 
the hypercube, and, in terms of absolute coordi- 
nates 

(11) 

with the overbar indicating forward averaging in 
the indicated coordinate, and _L indicating coor- 
dinates other than the one specified. The result 
of the above construction is that the link operator 
transforms very simply, 

= i 5 [<5^o(i M )i - (L^Sui], (12) 

where only the local value of the m M coordinate 
is shown (the other coordinates are the same 
throughout). The covariant, "transversely local" 
field strength is a path ordered product of link 
operators around the [i-v plaquette: 

exp(— igh 2 {T^ v ) m ) = P exp(— i<? <j> A ■ all) 

= (^) 10 (^)n(4)n(4)oi, (13) 

where now the local coordinates refer to the val- 
ues of m M and m v . This, however, does not re- 
duce to the free finite element form if g = 0, and 
therefore would violate unitarity. So we introduce 
a covariant averaging operator, 

(D X X) = ±1(L X ) 1 X 1 (L{) 1 +X \, (14) 

where now the coordinate refers to m\ . The field 
strength is then (symmetric averaging) 



(15) 



The field strength transforms covariantly in the 
sense 



(16) 



(£/Jm = exp(-i3/i(^) m ), 



(10) 



2.1. Yang-Mills equations 

Although the above field strength construction 
is local, it is not possible to obtain local Dirac 
and Yang-Mills equations without sacrificing uni- 
tarity H . The inductive approach in Q] gives the 
following covariant Yang-Mills equations 

E {{tt F ^ - ( F no}+^\=^, (17) 
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Table 1 



Ground-state energies for a Xq 2 k/2k potential. (Figures in parentheses are percentage relative errors.) 



k 




0(h) 




om 


Exact 


One state 


Two state 


One state 


Two state 


2 


0.4293(+2) 


0.4212(+0.1) 


0.4178(-0.7) 


0.4205(-0.06) 


0.42081 


3 


0.4639(+7) 


0.4391(+0.9) 


0.4453(+2) 


0.4328(-0.5) 


0.43493 


4 


0.5230(+13) 


0.4772(+3) 


0.5171(+11) 


0.4647(-0.2) 


0.46450 



where the value of the local v coordinate is shown 
(the curl involves averaging in mi), and the in- 
teraction term /C^ = K, n can be easily constructed 
from the difference equation 

2 

fcn-lL n + L n K n — J^L-n, {F tJ,l ')n,m ± }, (18) 

where L n — L v n . The current in ( |l7j ) may be 
taken to be the local form, 



3m = S^W^fe 
which transforms like (|lC 



(19) 



3. AXIAL- VECTOR ANOMALIES AND 
FERMION DOUBLING 

For an Abelian theory it is easy to obtain an 
explicit form for the transfer matrix, which carries 
the Dirac operator from one time to the next, 
VVi+i = T n ip n . In the temporal gauge, Aq = 0: 



T = 



1 - %i7°/2 + 7°7 • V 
l + «V7°/2-7°7-25' 



(20) 



where T> is a covariant derivative operator || . Us- 
ing this, nonzero vector and axial- vector anoma- 
lies have been computed in 2 and 4 space-time 
dimensions H|. For example, in 2 dimensions 



!"«"> 



e 2 E ( cos 2 tt/2M 
M sin7r/M \ 1 



(21) 



for M odd or even, respectively, where E is the 
electric field and M is the number of spatial lat- 
tice sites. Note this reduces to continuum form, 
with error, as expected, of 0(M~ 2 ), as M — * oo. 
At the same time we can compute the lattice vac- 
uum polarization (one fermion loop) and obtain 
results very close to the expected value of e 2 /7r. 
So there is no evidence of species doubling, which 



freedom is already signaled by the form of the free 
dispersion relation M, 



^ h 2 M 



(22) 



Note that the Minkowski-space formulation is 
crucial to this conclusion. If we were working in 
Euclidean space, it is possible to define an ac- 
tion from which the equations of motion could be 
derived. Since that action possess chiral symme- 
try, the currents derived from it (which involve all 
times, "past" and "future") will be anomaly free. 
Explicitly, fermion doublers then appear and con- 
tribute to the vacuum polarization |pj. 

But in Minkowski space, it is not possible to 
derive an action. Therefore, there is no connec- 
tion between symmetry and conserved currents. 
In fact, the current ([l9]) must be freely invented 
(consistent with gauge invariance), and can, and 
does, possess anomalies. Explicitly, the vacuum 
polarization gets its sole contribution from the 
neighborhood of zero-momentum fermions. 
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